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Most  Li  ion  insertion  batteries  consist  of  a  porous  cathode,  a  separator  filled  with  electrolyte  and  an 
anode,  which  very  often  also  has  some  porous  structure.  The  solid  part  especially  in  the  cathode  is  usually 
produced  by  mixing  a  powder  of  the  actual  active  particles,  in  which  Li  ions  will  be  intercalated,  binder 
and  carbon  black  to  enhance  the  electronic  conductivity  of  the  electrode.  As  a  result  the  porous  structure 
of  the  electrodes  is  very  complex,  leading  to  complex  potential,  ion  and  temperature  distributions  within 
the  electrodes.  The  intercalation  and  deintercalation  of  ions  cannot  be  expected  to  be  homogeneously 
distributed  over  the  electrode  due  to  the  different  transport  properties  of  electrolyte  and  active  particles 
in  the  electrode  and  the  complex  three-dimensional  pore  structure  of  the  electrode.  The  influence  of  the 
final  microstructure  on  the  distribution  of  temperature,  electric  potential  and  ions  within  the  electrodes 
is  not  known  in  detail,  but  may  influence  strongly  the  onset  of  degradation  mechanisms.  For  being  able 
to  numerically  simulate  the  transport  phenomena,  the  equations  and  interface  conditions  for  ion,  charge 
and  heat  transport  within  the  complex  structure  of  the  electrodes  and  through  the  electrolyte  filled 
separator  are  needed.  We  will  present  a  rigorous  derivation  of  these  equations  based  exclusively  on 
general  principles  of  nonequilibrium  thermodynamics.  The  theory  is  thermodynamically  consistent  i.e. 
it  guarantees  strictly  positive  entropy  production.  The  irreversible  and  reversible  sources  of  heat  are 
derived  within  the  theory.  Especially  the  various  contribution  to  the  Peltier  heat  due  to  the  intercalation 
of  ions  are  obtained  as  a  result  of  the  theory. 

Research  highlights:  ►  Thermodynamic  consistent  transport  theory  for  Li  ion  batteries  ►  Derivation  of 
all  irreversible  and  reversible  heat  sources  in  Li  ion  batteries  ►  Closed  set  of  equations  for  ion,  charge 
and  heat  transport  in  Li  ion  batteries  ►  Theory  of  Peltier  heat  for  Li  ion  intercalation  ►  Microstructure 
resolved  transport  in  porous  electrodes 
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1.  Introduction 

Understanding  the  heat  generation  and  temperature  distribu¬ 
tion  in  Li  ion  battery  cells  is  of  great  practical  interest.  Many 
degradation  processes  in  Li  ion  batteries  are  caused  by  and  most 
of  them  are  enhanced  under  temperature  increase  [1,2].  Thermal 
runaway  may  begin  as  a  very  local  phenomenon  by  initiating  an 
exothermal  reaction  at  a  local  hotspot.  Therefore  the  average  tem¬ 
perature  of  a  cell  is  not  sufficient  to  determine  the  beginning  of 
a  thermal  runaway.  A  spatially  resolved  mathematical  description 
of  the  processes  on  the  scale  of  the  cell  would  allow  to  simulate 
the  behavior  of  the  cell  under  varying  operating  conditions  and  to 
detect  the  probability  of  dangerous  hotspots. 

Mathematical  modeling  of  Li-ion  batteries  on  cell  resolved  level 
was  pioneered  by  the  work  of  Newman  and  his  coworkers  [3-5] 
and  extended  and  refined  by  many  other  authors  [6-8].  The  mod¬ 
eling  approaches  are  based  on  transport  equations  for  Li  ions  and 


*  Corresponding  author.  Tel.:  +49  631  316004301;  fax:  +49  631  316005301. 
E-mail  address:  Arnulf.Latz@itwm.fraunhofer.de  (A.  Latz). 

0378-7753 /$  -  see  front  matter  ©  2010  Elsevier  B.V.  All  rights  reserved, 
doi:  10.1 01 6/j.jpowsour.201 0.1 1.088 


charges  in  the  electrolyte  as  well  as  in  the  active  particles  of  cathode 
and  anode  under  isothermal  conditions.  Most  modeling  approaches 
for  the  thermal  behavior  of  batteries  are  concentrating  on  overall 
thermal  balance  equations  for  a  whole  cell  [9-14]  by  combining 
phenomenologically  thermodynamic  considerations  on  entropy  or 
enthalpy  changes  within  a  cell  with  reasonable  assumptions  on 
out  of  equilibrium  processes  like  Joule  heating,  heat  of  mixing, 
Peltier  effect  and  Soret  effect.  An  approach  based  on  local  continu¬ 
ity  equations  for  the  temperature  was  presented  in  [  1 5],  extending 
previous  work  on  species  and  charge  transport  in  batteries  [16]. 
These  authors  focused  on  deriving  macroscopic  equations  for  the 
heat  transport  in  porous  electrodes  using  the  volume  averaging 
technique.  The  microscopic  equation  for  the  temperature  was  for¬ 
mulated  without  considering  possible  changes  of  the  transport 
equations  describing  species  and  charge  transport.  Also  the  inter¬ 
face  conditions  have  not  been  derived.  Fully  coupled  models  have 
been  considered  for  battery  stacks  [17]  and  other  types  of  batter¬ 
ies  (see  e.g.  [18,19]).  Also  in  these  cases  the  equations  were  not 
derived  but  formulated  as  balance  equations  in  which  known  phys¬ 
ical  effects  as  well  as  the  reaction  kinetics  for  all  relevant  chemical 
reactions  were  phenomenologically  incorporated. 
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A  different  approach  to  treat  transport  in  electrochemical 
systems  is  the  systematic  approach  of  nonequilibrium  thermo¬ 
dynamics  [20]  pioneered  by  Landau  [21,22].  In  this  approach 
only  conservation  equations  and  general  principles  of  statistical 
mechanics  were  used  to  derive  the  general  form  of  the  relevant 
transport  equations.  Only  measurable  quantities  like  thermody¬ 
namic  derivatives  and  transport  coefficients  enter  the  equations. 
The  Onsager  relations  which  guarantee  the  positivity  of  the  entropy 
production  and  impose  strict  relations  between  the  equations  for 
densities,  charges  and  temperature,  play  a  central  part  in  the 
derivation.  Chemical  reactions  can  easily  be  treated  within  this 
formalism  [20].  The  interface  conditions  can  also  be  derived  by 
necessary  continuity  constraints  (as  e.g.  conservation  of  charges, 
energy,  etc.)  and  considering  the  physical  and  chemical  pro¬ 
cesses  within  the  interface  [20].  Contrary  to  the  phenomenological 
approaches  the  terms  contributing  to  the  heat  production  are 
systematically  derived  using  nothing  else  but  generally  valid  ther¬ 
modynamic  relations  and  total  energy  conservation,  instead  of  just 
physical  intuition  and  experimental  experience.  In  this  manner  it 
is  made  certain  that  no  source  of  heat  is  neglected,  if  all  hydrody¬ 
namic  variables  and  all  relevant  chemical  reactions  are  taken  into 
account. 

An  especially  elegant  form  of  this  formalism  was  recently  used 
to  derive  a  hydrodynamic  theory  of  electromagnetic  field  in  con¬ 
tinuous  media  [23-26],  which  can  also  be  applied  to  charged 
and  magnetic  liquids  [27].  Also  charge  and  species  transport  in 
Li  ion  batteries  can  be  dealt  with  using  the  nonequilibrium  ther¬ 
modynamics  framework  [28].  We  will  derive  in  the  following  the 
complete  set  of  equations  including  all  interface  conditions  for  the 
ion,  charge  and  thermal  transport  in  Li  ion  battery  cells.  We  will 
consider  a  spatial  scale  L  which  is  large  compared  to  the  scale  Ldl  of 
the  diffuse  double  layer  (i.e.  ~10-20nm  [4]).  Concentrated 

electrolytes  are  neutral  on  this  scale  due  to  the  strong  Coulomb 
interactions  between  the  ions  leading  to  an  extremely  small  Debye 
screening  length.  Only  in  the  double  layer  around  the  active  parti¬ 
cle,  which  usually  have  diameters  on  micrometer  scale,  potential 
gradients  are  strong  enough  to  cause  charge  separation  [4].  Also 
active  particles  fulfill  charge  neutrality  although  the  mechanism 
for  obtaining  charge  neutrality  is  very  different.  In  active  particles 
it  is  the  large  mobility  of  electrons  which  guarantees  charge  neu¬ 
trality.  The  charge  of  an  inserted  Li  ion  is  instantaneously  shielded 
by  local  rearrangements  of  electronic  charges  and  the  transport  of 
electrons  into  the  active  particles  over  the  current  collectors.  In  the 
electrolyte  charge  neutrality  will  lead  to  highly  correlated  motion 
of  negative  and  positive  ions.  We  will  make  explicit  use  of  charge 
neutrality  in  the  derivation  of  our  equations. 

In  addition  to  naturally  reducing  the  numbers  of  relevant 
equations,  the  use  of  charge  neutrality  also  guarantees,  that  the 
transport  coefficients  appearing  in  theory  are  the  measurable  quan¬ 
tities  of  the  electrolyte.  Let  us  for  example  consider  the  diffusion 
processes  in  a  molten  salt  consisting  of  positive  and  negative  ions. 
In  general  the  self-diffusion  coefficients  are  different  for  the  two 
kinds  of  ions.  If  there  were  no  interaction  between  the  ions,  we 
could  expect  to  obtain  free  electric  charges  in  relaxation  processes 
even  if  the  initial  perturbation  were  neutral  due  to  the  different 
diffusion  length  of  the  ions.  In  reality  the  strong  Coulomb  interac¬ 
tion  between  the  ions  prevents  the  appearance  of  free  charges.  The 
process  relevant  for  diffusion  is  in  fact  the  collective  interdiffusion 
process  with  a  uniquely  defined  interdiffusion  coefficient  for  both 
species.  This  unique  quantity  for  positive  and  negative  ions  will 
appear  naturally  in  our  theory. 

To  derive  the  transport  equations  in  the  electrolyte  and  the 
solid  particles  we  use  a  general  and  rigorous  theory  for  polarizable, 
conducting  media  introduced  in  [23]  and  generalized  in  [24].  The 
theory  is  based  on  general  thermodynamic  principles,  Maxwells 
theory  for  electromagnetic  fields  and  Onsagers  reciprocal  relations. 


Although  transport  mechanisms  in  active  particles  and  electrolyte 
can  be  very  different  on  the  microscopic  (atomic)  scale,  their  macro¬ 
scopic  form  just  differ  in  the  functional  dependencies  and  the  size 
of  the  transport  coefficients.  For  example  in  active  particles  charge 
is  transported  mainly  by  pure  electronic  conduction.  The  contri¬ 
bution  of  the  ion  diffusion  in  the  active  particles  to  the  electric 
current  can  be  neglected  due  to  the  large  mobility  of  the  electrons 
compared  to  the  ions.  Charge  transport  in  the  electrolyte  on  the 
other  side  is  exclusively  due  to  ionic  transport.  In  fact  the  transfer 
of  electrons  into  the  electrolyte  would  result  in  the  reduction  of  Li 
ions  in  the  electrolyte  to  metallic  lithium  and  is  considered  to  be 
one  of  the  many  degradation  mechanisms  in  Li  ion  batteries  [1]. 
Macroscopically,  the  differences  in  the  microscopic  charge  trans¬ 
port  mechanisms  are  reflected  in  different  functional  dependencies 
of  the  charge  conductivity  coefficients  on  temperature  and  Li  ion 
concentration.  The  ionic  conductivity  in  the  electrolyte  is  usually 
strongly  dependent  on  the  concentration  of  Li  ions.  The  electronic 
conductivity  of  the  active  particles  is  basically  not  affected  by  it  and 
can  be  assumed  to  constant.  In  addition,  the  transference  number 
of  ions  in  active  particles  is  approximately  zero  due  to  the  vanishing 
contribution  of  ion  diffusion  to  the  electrical  current. 

2.  Model 

Originally  the  theory  of  Henjes  and  Liu  [23]  was  used  to  derive 
the  hydrodynamic  equations  for  a  one  component  polarizable  liq¬ 
uids.  But  the  principles  can  easily  be  applied  to  a  mixture  of  a 
dissociating  salt  and  a  solvent  i.e.  for  a  generic  electrolyte  as  used 
in  batteries.  The  reduction  to  the  transport  in  the  solid  particle  is 
straight  forward  after  the  full  set  of  equations  for  the  electrolyte 
are  derived.  Although  we  have  to  deal  with  a  at  least  three  compo¬ 
nent  system  of  positive  and  negative  ions  as  well  as  neutral  solvent 
molecules,  the  problem  can  be  reduced  to  an  effective  one  com¬ 
ponent  problem  using  two  constraints,  which  are  applicable  to 
the  physical  situation  of  a  concentrated  electrolyte  in  a  battery. 
The  first  observation  is  that  convective  transport  is  highly  unlikely 
under  normal  operating  conditions  in  a  battery.  We  therefore  may 
safely  assume  that  the  center  of  mass  of  any  hydrodynamic  vol¬ 
ume  element  is  at  rest.  This  assumption  allows  to  eliminate  the 
concentration  of  the  neutral  solvent  as  independent  variable.  With 
M0,  M+,  M_  being  the  molar  masses  of  solvent  and  positive  and 
negative  ions  respectively  we  get  in  the  absence  of  convection  the 
relations: 

M0dc0  +  M+dc+  +  M_dc_  =0  (1 ) 

for  the  changes  in  the  respective  concentrations. 

The  second  observation  concerns  the  occurrence  of  finite  free 
charges  in  a  volume  element.  It  can  easily  be  estimated  [4]  that 
the  electric  fields  which  where  necessary  to  impose  charge  sep¬ 
aration  on  a  micrometer  scale  are  much  bigger  than  the  average 
fields  present  in  a  battery.  Charge  separation  is  only  observed  in 
the  double  layer  around  active  particles  on  a  scale  of  the  order 
of  10-20  nm.  The  present  theory  will  be  restricted  to  scales  above 
100  nm.  A  detailed  theory  for  the  linear  and  nonlinear  phenomena 
in  the  double  layer  can  be  found  in  [29,30].  If  we  denote  molar  con¬ 
centrations  of  positive  and  negative  ions  of  charge  z+  and  z_  with 
c+  and  c_,  respectively,  charge  neutrality  requires 

v+z+c+  +  V-Z-C _  =  0  (2) 

where  v+,  v_  are  the  stoichiometric  numbers  of  positive  and 
negative  charges.  It  is  therefore  sufficient  to  determine  the  trans¬ 
port  equations  for  the  concentration  c  =  c+ =  -  (v_z_/v+z+)c_.  For 
notational  simplicity  we  restrict  our  theory  to  a  binary  salt  with 
v+  =  V-  =  \  and  consequently  z_  =  -  z+  as  well  as 

C+  =  C_ . 


(3) 
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The  fundamental  thermodynamic  relation  for  a  polarizable  system  difficult  to  measure  it  is  more  convenient  to  work  with  c,  T  and 

in  an  external  electromagnetic  field  is  given  by  [23]  E  as  independent  variable.  We  therefore  rewrite  Eq.  (16)  as 


de  =  Tds  +  fi+dc+  +  /x_dc_  +  /x0dc0  +  E  ■  dD  +  H  dB  (4) 

Here  e,  s ,  ii[  are  the  energy  density,  entropy  density  and  the 
chemical  potential  of  species  i.  E,  D,  B,  H  are  electric  field,  electric 
displacement,  magnetic  induction  and  magnetic  field,  respectively. 
Using  (1 )  and  (3),  Eq.  (4)  reduces  to 

de  =  Tds  +  fide  +  £  •  db  +  H  •  dB  (5) 

where  c  =  c+  is  the  concentration  of  Li-ions  and  the  effective  chem¬ 
ical  potential  /n  of  the  lithium  ions  is  given  by 


R  =  -Q  ^~  (^jN+  Vc+j-E  (17) 

with  Q  =  q  +  N+T(3/x/9T).  The  expression  for  the  entropy 
production  defines  ^  the  independent  thermodynamic  fluxes 
Q/T,  (9/x/9c)N+  and  j,  for  which  constitutive  equations  have  to  be 
formulated.  In  order  to  guarantee  positiveness  of  the  entropy  pro¬ 
duction  R ,  the  constitutive  equations  have  to  be  chosen  such,  that 
the  expression  for  R  ( 1 7)  is  a  positive  definite  binary  form.  The  most 
general  expression  for  the  fluxes  is  therefore  given  by 


M+  z+  (  M_  \ 

For  a  binary  charge  neutral  salt  this  reduces  to 

M_  +  M+ 

=  +  - — — fi0 


(6) 

(7) 


The  next  building  blocks  are  the  conservation  equation  for  con¬ 
centration  of  Li  ion,  charge  density  pd ,  entropy  as  well  as  the 
Maxwell  equations.  Note  that  the  macroscopic  velocity  is  set  to 
zero,  due  to  the  main  assumption  of  convection  free  transport  in  Li 
ion  batteries.  The  equation  for  pd  is  obtained  from  the  conservation 


equations  of  ions: 

t> 

i 

II 

(8) 

dc  ->  -> 

(9) 

—  =-V  N- 

Here  N+,  N_  are  the  fluxes  (i.e.  the  number  of  moles  per  unit  time 
and  unit  area)  of  the  positive  Li  ions  and  the  negative  counter  ions, 
respectively.  Using  (3)  we  obtain  with  pd  =  z+c+  +z_c_  and  the  def¬ 
inition  of  the  electric  current  j  =  z+N+  +  z_N_ : 


-^r  =  0  =  -  V  -  j  (10) 

From  the  Maxwell  equation  follows: 

E  ■  dtD  +  H  ■  dtB  =  -V  ■  (E  x  H)  -  j  E  (11) 


The  general  form  of  the  entropy  equation  is  given  by 


Here  q  is  the  heat  flux.  The  entropy  production  R  can  be  derived  by 
combining  the  first  law  (5)  and  the  conservation  of  the  total  energy. 
From  the  first  law  (5)  we  obtain: 


dte  =  Tdts  +  fjidtc  +  EdtD  +  HdtB  (13) 

With  (8),  (11)  and  (12),  Eq.  (13)  transforms  into 

-  _  VI  -  -  - 

dte  =  -V.(q  +  N+/x  +  FxH)  +  K  +  q.  —  -j-E  +  N+-  V/z  (14) 

Due  to  the  conservation  of  total  energy  the  change  in  the  local 
energy  has  to  fulfill  a  continuity  equation  i.e.: 

3te  =  -V  Je  (15) 


By  comparing  (14)  and  (15)  we  obtain  the  entropy  production  R  as 

R  =  -q.fl-N+.Vn+j.E  (16) 

The  entropy  production  has  to  be  strictly  positive  or  zero.  This 
requirement  imposes  strong  constraints  on  the  constitutive  rela¬ 
tions  for  the  fluxes  q ,  N+  and  j.  Since  the  chemical  potential  is 


-T  7 


dfl  \  Kt+ 
~dc)Ez~+ 


-/k 


-Pk 

y 

T 


(18) 


The  symmetry  of  the  Onsager  matrix  in  Eq.  (18)  is  necessary  for 
the  positiveness  of  the  entropy  production  and  a  consequence  of 
fundamental  principles  of  nature  [20].  The  physical  meaning  of  the 
various  transport  coefficients  a ,  /c,  y,  ft,  v,  t+  is  getting  more  obvious 
after  transforming  equivalently  the  constitutive  equations  into  the 
more  traditional  form: 


N+  =  -DV c  +  ~pl  ~  VT 

a— 

In  this  form  the  entropy  production  has  the  simple  form: 


T  k 


( N+-(t+/Fz+)j ) 
dc 


D 


(19) 


(20) 


(21) 


(22) 


In  order  to  be  positive  each  term  in  the  expression  (22)  has  to  be 
positive  i.e.  the  heat  conductivity  k,  the  electric  conductivity  k  and 
the  interdiffusion  coefficient  D  have  to  be  positive. 

Eq.  (20)  can  also  be  expressed  using  the  electrochemical  poten¬ 
tial  (p  =  fil(z+F)  +  O,  where  O  is  defined  by  E  =  -  VO: 


j  =  -KV(p  -  K 


4-1 

z+F 


1 


z+F  \  8T 


^  I  I kVT 


(23) 


The  electric  conductivity  k  in  the  electrolyte  is  a  pure  ionic  conduc¬ 
tivity.  In  addition  to  the  ionic  conductivity,  the  transference  number 
t+  of  Li  ions  and  the  Seebeck  coefficient  are  the  relevant  transport 
coefficients  for  the  electric  current.  The  transference  number  also 
appears  in  the  relation  for  the  ionic  flux  N+  and  determines  the 
amount  of  ion  flux  contributed  by  the  electric  current.  The  Seebeck 
coefficient  has  the  units  VK-1  and  is  related  to  the  Peltier  coef¬ 
ficient  n  by  Tl  =  T/3.  It  quantifies  the  thermoelectric  effect  in  the 
electrolyte  or  the  active  particles.  The  Seebeck  coefficient  can  in 
principle  be  positive  or  negative. 

The  two  additional  transport  coefficients  relevant  for  the  ionic 
flux  are  the  inter  diffusion  coefficient  D  and  the  Soret  coefficient 
kj.  The  interdiffusion  coefficient  D,  appearing  in  the  relation  for  the 
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ionic  flux  N+,  is  measured  in  the  limit  of  zero  electrical  current  and 
zero  temperature  gradient.  It  is  given  by 


Since  D  has  to  be  positive  the  Onsager  coefficient  a  and  the  trans¬ 
ference  number  fulfill  an  obvious  inequality  relation.  The  Soret 
coefficient  kT  quantifies  thermodiffusional  processes  and  is  given 
by 


D  c  kj 
T 


c p± 

Fz+ 


(25) 


Usually  the  Soret  coefficient  is  of  order  0.6-1 .5  [31  ],  but  the  over¬ 
all  size  of  the  Soret  effect  is  believed  to  be  very  small  in  liquids 
and  electrolytes  [32].  Therefore  the  ion  flux  induced  by  tempera¬ 
ture  gradients  is  neglected  compared  to  the  one  caused  by  electric 
currents  and  gradients  in  concentrations. 

The  last  transport  coefficient  is  the  thermal  conductivity  X  mea¬ 
sured  in  the  limit  of  vanishing  electric  current  and  vanishing  ionic 
flux.  It  is  given  by 


A  =  K“(t^)  f^c)2-^2 


(26) 


It  also  has  to  be  strictly  positive  in  order  to  guarantee  positive 
entropy  production  R. 


2.1.  Heat  transport 


The  temperature  changes  due  to  thermal  conduction  and  four  dif¬ 
ferent  sources  for  heat.  They  are  in  the  order  of  their  appearance 
in  Eq.  (31):  Joule’s  heat,  Thompson  effect,  heat  of  mixing  and  the 
Soret  effect.  The  Thompson  effect  can  also  be  written  as 

TV-($)  =  7-My.VT  :=Mr]  vr  (32) 


Here  we  used  charge  conservation  (Eq  (10))  and  the  Thompson  rela¬ 
tion  for  the  Thompson  coefficient  /jlt  =  T{  d P/dT).  If  we  neglect  all 
contributions  proportional  to  kT  the  heat  Eq.  (31)  reduces  further 
to 


cppdtT  =  V  •  (AVI)  + 


l 

K 


Tv-m+ 


(33) 


The  equations  for  ion  concentration,  electric  field  respective  electric 
potential  defined  by  E  =  -VO  and  temperature  can  be  applied  in 
the  electrolyte  as  well  as  in  the  active  particles.  The  only  difference 
is  the  electric  conductivity,  which  is  replaced  by  the  average  elec¬ 
tronic  conductivity  o  of  the  solid  active  particle  and  carbon  black, 
which  is  usually  used  to  enhance  the  electronic  conductivity  of  the 
electrodes. 


2.2.  Interface  conditions 

To  couple  the  transport  in  electrolyte  and  active  particles  we 
have  to  formulate  interface  conditions,  which  describe  the  interca¬ 
lation  and  deintercalation  of  ions. 


The  equation  for  the  temperature  can  be  derived  from  the 
entropy  balance  (12)  and  the  expression  for  the  entropy  production 
(17)  and  the  relation: 

TdtS  =  CppdtT +  T^\Tdtc  (27) 

where  cp  and  p  are  the  specific  heat  per  unit  mass  and  the  mass 
density  of  the  electrolyte  or  the  active  particles.  Using  the  thermo¬ 
dynamic  relation (  9s/9c)  |  T  =  -{d /i/dT)  |  cand  the  continuity  Eq.(8) 
we  obtain: 


cppdj  =  -V-q-N+-V/i+j  E-T(  ^  )  V  ■  N. 


(28) 


Transforming  from  the  original  heat  flux  q  to  the  renormalized  heat 
flux  Q.  we  get 


CppdtT  —  — V  •  (  Q.  -  TN+  (  ^ 


-iSKw-’ft 


V-N, 


(29) 


The  third  and  last  term  on  the  right  hand  side  of  Eq.  (29)  are  can¬ 
celled  by  the  first  one  on  the  right  hand  side,  if  we  neglect  the 
spatial  variation  of  the  thermodynamic  derivative  (d/u/dT).  The 
temperature  equation  is  therefore  given  by 


CppdtT  =  -V  •  Q  -  (  ^ )N+Vc+j-E 


(30) 


Using  the  constitutive  relations  (19)-(21),  Eq.  (30)  can  be  trans¬ 
formed  into 


cppdtT  =  V.(AVT)  +  J--TV.(^)+  ( 


dfi\  ( N+-(t+/Fz+)j ) 


2.2.1.  Interface  conditions  for  ionic  flux  and  electric  current 

The  interface  conditions  describe  the  intercalation  reaction  and 
the  deintercalation  reaction  respectively  on  the  mesoscopic  scale 
(i.e.  beyond  the  scale  of  the  diffuse  layer  [4]).  For  one  step  reactions 
it  is  usually  assumed  that  the  transport  of  ions  across  the  interface 
is  completely  described  by  the  Butler  Volmer  approach.  More  com¬ 
plicated  reactions  may  require  the  use  of  more  elaborate  theories 
[33],  which  can  be  easily  incorporated  within  our  theory.  Here  we 
use  for  simplicity  the  Butler  Volmer  theory.  The  current  density 
across  the  interface  ise  due  to  the  intercalation  reaction  is  within 
the  Butler  Volmer  approach  given  by  [4]: 


io  (e 


ise  =  i0  exp 


(XaF 

RT 


As 


-  exp 


-acF 

~ RT 


ris  ) 


(34) 


aA  and  ac  with  aA  +  ac  =  1  are  weighting  the  anodic  and  the  cathodic 
contribution  of  the  overpotential  rjs  to  the  overall  reaction.  The 
overpotential  is  the  deviation  of  the  electrochemical  potential  from 
the  chemical  equilibrium  between  active  particle  and  electrolyte. 
It  is  therefore  defined  by 


J)S  :=  Os  -  Oe  -  (35) 

The  overpotential  vanishes  obviously  if  the  active  particle  is  in  equi¬ 
librium  with  the  electrolyte.  Usually  the  overpotential  is  expressed 
with  the  help  of  the  half  cell  open  circuit  potential  U0  of  the 
respective  electrode  relative  to  a  Li  metal  electrode.  Without  loss 
of  generality  setting  0Ll  =  0,  the  open  circuit  potential  U0  can  be 
written  as 

U0  =  (36) 

Z+r 

Using  this  expression  in  Eq.  (35)  we  get 


ijs  :=  Os -  Oe -  U o  - 


Me  -  Mli 


z+F 


(37) 


-TV- 


Since  the  chemical  potential  of  the  electrolyte  /xe  is  very  different 
from  the  chemical  potential  /xLl  of  the  Li  metal  electrode  the  last 
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two  terms  in  Eq.  (37)  do  not  cancel  in  general.  Introducing  the  elec¬ 
trochemical  potential  of  the  electrolyte  cpe  relative  to  the  chemical 
potential  of  Li  -  metal 

(Pe  :=  <t>e  +  (38) 

z+t 

we  can  rewrite  (37)  as 

Vs  :=  -  <Pe  -  U o  (39) 


left  hand  side  of  Eq.  (46)  vanishes.  Using  Eqs.  (36)  and  (38)  we 
obtain: 

ft  ■  (4  -  Qe)  =  -isetls  ~  iseT ^  (47) 

Neglecting  the  temperature  dependence  of  the  chemical  potential 
of  Li  metal  and  using  the  constitutive  relation  for  Q  Eq.  (21)  we 
finally  obtain: 


and  obtain  the  standard  form  of  the  overpotential  but  with  the 
electrochemical  potential  instead  of  the  electric  potential. 

The  amplitude  i'o  in  Eq.  (34)  is  given  by 

i0  =  kc““c“o(cs,max-cs)“c  (40) 


k  is  a  reaction  rate.  cs,max  is  the  maximum  concentration  of  ions  in 
the  active  particle.  We  assume  that  Li  ions  are  not  stored  in  the 
double  layer  (i.e.  all  Li  ions  are  intercalated  in  the  active  particle 
or  released  into  the  electrolyte).  There  should  also  be  no  flux  of 
negative  charges  across  the  double  layer.  Neither  enter  electrons 
the  electrolyte  nor  intercalate  negative  ions  from  the  electrolyte  in 
the  active  particles  under  ideal  conditions.  Both  effects  would  lead 
to  degradation  [1  ],  which  could  of  course  be  modeled  by  modifying 
the  interface  conditions.  The  absence  of  flux  of  negative  charges 
across  the  interface  especially  means  that  the  total  current  across 
the  electrolyte-particle  interface  is  due  to  transport  of  positive  ions 
only.  If  the  particle  is  completely  filled  i.e.  c  =  cs,max ,  it  has  to  be  made 
sure  by  the  interface  conditions  that  no  electrical  current  j  is  carried 
by  negative  charge  carriers  across  the  interface.  These  conditions 
can  be  formulated  mathematically  in  the  following  way  with  the 
normal  h  pointing  from  the  solid  into  the  electrolyte: 


—  he 

(41) 

=  he 

(42) 

z+F 

(43) 

z+F 

(44) 

-Asn  •  VTS  +  Xeh  ■  VTe 


+  ise 


dU0. 

CS1FI<T,S+Ce 


dfte  fcr,e(l  -  4)\ 

dc  Z+F  J 


(48) 


Here  we  also  used  that  the  transference  number  of  Li  ions  in  the 
active  particle  may  be  neglected  due  to  the  high  mobility  of  elec¬ 
trons  compared  to  the  mobility  of  Li  ions.  The  physical  meaning  of 
the  various  expressions  on  the  right  hand  side  of  Eq.  (48)  is  obvi¬ 
ous.  The  first  one  is  the  irreversible  heat  production  due  to  Joule 
heating.  The  next  term  contributes  to  the  reversible  Peltier  effect 
and  the  last  one  is  the  Soret  effect  in  the  double  layer.  The  Peltier 
coefficient  is  given  by 


rr  t(  a  r  \  i  +  /W(z+F)) 

U  =  T(PS  -  Pe)  +  T - — - 


(49) 


Note  that  in  [11]  only  the  partial  derivative  of  the  open  circuit 
potential  was  given  as  Peltier  coefficient.  Our  rigorous  approach 
shows  naturally  that  also  the  differences  in  the  Seebeck  coefficients 
of  the  two  phases  and  the  thermal  derivative  of  the  electrolytes 
chemical  potential  are  contributing  to  the  Peltier  coefficient.  Since 
our  approach  is  based  on  experimentally  accessible  transport 
coefficients  it  is  either  possible  to  measure  directly  the  Peltier  coef¬ 
ficient  [34]  or  to  deduce  it  from  measurements  of  the  open  circuit 
potential,  the  Seebeck  coefficients  of  the  two  phases  [35]  and  the 
thermodynamic  derivatives  of  the  chemical  potential  of  the  elec¬ 
trolyte. 


2.2.2.  Thermal  interface  conditions 

To  derive  the  thermal  interface  conditions  we  consider  a 
infinitesimally  extended  small  piece  of  the  interface  perpendicu¬ 
lar  to  the  normal  h  which  points  from  the  active  particle  into  the 
electrolyte.  The  interface  conditions  are  most  conveniently  derived 
using  the  form  (28)  of  the  heat  transport  equation.  We  integrate  the 
temperature  balance  equation  over  the  infinitesimal  small  volume 
element,  which  contains  the  whole  thickness  of  double  layer.  The 
term  proportional  to  V  •  N+  does  not  contribute  to  the  flux  since  the 
flux  is  conserved  across  the  double  layer  (see  Eq.  (43)).  We  obtain 


-// 


dVcppdtT  =  /  /  dA(n  ■  (qs  -  qe)  +  n  ■  N+(ns  -  p,e) 


+  n -J(<J>S  -  <I>e)) 

Using  Eqs.  (41 )  and  (43)  and  the  definition  of  Q  we  obtain 

JJJ  dVcppdtT  =  jj  dA(ii  ■  (4  -  Qe)) 

+  -  /te  +Z+F(< I>S  -  <f>e)) 

Z+t 


(45) 


he  j 

z+F 


9/Xs 

dT 


9/Xe 

dr 


(46) 


If  we  assume  that  the  double  layer  instantaneously  changes 
temperature  upon  changes  of  the  currents,  concentrations  and  tem¬ 
peratures,  the  double  layers  is  always  in  a  stationary  state  i.e.  the 


2.3.  Closed  set  of  equations  for  transport  in  batteries 


2.3.1.  Electrolyte 

The  natural  variables  for  the  electrolyte  are  the  ion  concentra¬ 
tion  c,  the  electrical  potential  and  the  temperature  T.  Due  to 
the  experimentally  motivated  formulation  of  the  interface  con¬ 
ditions  with  the  open  circuit  potential  U0  relative  to  a  Li  metal 
electrode  instead  of  the  difference  between  chemical  potentials 
of  solid  particle  and  electrolyte  the  electrochemical  potential  <p 
(see  Eq.  (38))  had  to  be  introduced.  It  is  therefore  more  conve¬ 
nient  to  formulate  the  transport  equations  with  <pe  insteadt  of  Oe. 
One  has  to  keep  in  mind  that  the  electric  field  is  still  given  by 
E  =  -VOe  =  -V(pe  +  V/xe/(z+F).  The  transport  equations  are  then 
given  by 

3 tCe  =  v  •  (DeVCe)  -  V  •  (3±j)  +  V  .  (^Ivt)  (50) 

0  -V 

+v.Wft-jE(|j)vn 


Cp.eMT  =  V  •  (A.eVT)  +- - TV  • 


m-£j> 


-TV- 


HD 


(52) 
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Note  that  the  first  two  equations  without  the  terms  proportional 
to  the  temperature  gradients  have  the  form  of  the  equations  for 
the  concentrated  electrolyte  by  Newman  [4],  but  it  is  important 
to  realize  that  there  exists  an  important  quantitative  difference. 
In  the  concentrated  electrolyte  theory  of  [4]  the  potential  in  Eq. 
(51)  is  interpreted  as  the  electrical  potential  Oe.  In  our  theory 
the  electrical  potential  had  to  be  replaced  by  the  electrochemical 
potential  <pe  of  the  electrolyte,  if  we  use  the  open  circuit  potential 
relative  to  Li  metal  in  the  simple  expression  for  the  overpotential 
Eq.  (39). 

2.3.2.  Active  particles 

Our  derivation  of  the  constitutive  relations  for  concentration, 
electrical  potential  and  temperature  is  based  on  general  thermo¬ 
dynamic  principles  and  can  therefore  be  applied  equally  to  the 
transport  in  active  particles.  Only  if  the  transport  is  anisotropical 
the  transport  coefficients  have  to  be  replaced  by  tensorial  quan¬ 
tities.  Here  we  restrict  ourselves  to  active  particles  which  behave 
approximately  isotropic  on  the  micrometer  scale.  The  equations 
are  slightly  simplified  due  to  the  reasonable  assumption  that  the 
transference  number  can  be  set  to  zero,  since  the  electrical  current 
is  predominantly  electronic  of  nature.  We  will  denote  the  electronic 
conductivity  with  a.  The  transport  equations  are  then  given  by 

dtCs  =  V  ■  ( DsVcs )  +  v  ■  ^DCsrkr’s  vrj  (53) 

0  =  V-(crV0s)  +  V-(/Ss<rVT)  (54) 

cP,spdtT  =  V  •  (X,VT)  +  —  TV  •  (M  -  z+F^° 

+  Tz+FV(cs^.'^N+)  (55) 

Together  with  the  interface  conditions  (41 )-  (44)  and  (48)  the  com¬ 
plete  set  of  equations  for  transport  in  Li-ion  batteries  is  formulated. 
Of  course  the  equations  have  to  be  closed  by  using  appropriate 
boundary  conditions,  but  this  is  standard  procedure.  They  may 
change  from  application  to  application  depending  on  how  the  cell  is 
coupled  to  the  outside  world.  Especially  different  cooling  procedure 
will  require  different  boundary  conditions  either  for  the  thermal 
flux  or  the  temperature  itself.  Also  the  coupling  of  the  current 
collector  to  the  electrodes  will  influence  the  choice  of  boundary 
conditions. 

3.  Conclusions 

We  have  derived  the  complete  set  of  transport  equations  on  a 
spatial  scale  larger  than  the  diffuse  double  layer  in  Li  ion  batteries. 
Our  derivation  is  based  on  general  principles  of  nonequilibrium 
thermodynamics,  which  guarantee  thermodynamic  consistency 
and  especially  strictly  positive  entropy  production.  The  equa¬ 
tions  for  transport  in  the  electrolyte  differ  from  the  concentrated 
electrolyte  theory  of  [4].  It  can  easily  be  shown  [28],  that  the  con¬ 
centrated  electrolyte  theory  of  [4]  violates  the  Onsager  relations. 
But  we  have  also  shown  above,  that  the  coupled  equations  for 
charge  and  species  transport  in  Li  ion  batteries  (neglecting  heat 
transport)  can  be  transformed  in  the  form  used  in  [3],  [4],  if  in  the 
electrolyte  the  pure  electrical  potential  is  replaced  by  the  electro¬ 
chemical  potential  in  combination  with  the  definition  (39)  for  the 
overpotential. 

The  complete  set  of  equations  can  either  be  used  to  simulate 
the  transport  within  cells  using  a  spatial  representation  of  the  elec¬ 
trodes,  which  resolves  the  microstructure  of  the  porous  electrodes 
or  as  a  starting  point  to  obtain  a  porous  media  representation  of 
the  electrodes  as  in  [15].  In  the  microscopic  theory  active  particles 


and  electrolyte  are  treated  as  separate  media.  So  far  this  approach  is 
only  used  for  simulating  the  ion  and  charge  transport  under  isother¬ 
mal  conditions  [36,37].  Inclusion  of  heat  transport  is  planned  as 
future  work.  The  computational  complexity  of  the  problem  allows 
only  to  simulate  representative  volume  elements  (REVs)  of  bat¬ 
tery  cells.  But  as  post  processing  effective  properties  of  the  porous 
electrode  may  be  obtained  as  e.g.  effective  diffusion  coefficients  or 
conductivities  by  numerically  averaging  over  the  REV.  These  quan¬ 
tities  may  then  be  used  for  simulating  the  whole  cell  modeled  with 
the  porous  electrode  model. 

So  far  side  reactions  as  e.g.  formation  of  the  solid-electrolyte 
interface  (SEI)  or  transport  within  a  SEI  have  not  been  considered. 
Within  our  approach  it  is  straight  forward  to  include  those  phenom¬ 
ena.  Surface  reactions  will  lead  to  a  modification  of  the  interface 
conditions.  Transport  within  a  SEI  requires  the  addition  of  another 
thermodynamic  phase  with  appropriate  transport  coefficients  [38]. 
Volumetric  reactions  require  the  inclusion  of  the  reacting  species 
and  the  information  on  the  reaction  kinetics  [20].  The  derivation 
of  the  transport  equations  will  proceed  along  the  same  line  as 
shown  above  with  the  only  difference  that  the  concentrations  of 
the  species  are  not  conserved  separately  but  just  the  sum  of  them. 
In  future  work  we  will  make  extensive  use  of  the  model  to  obtain 
a  better  understanding  of  the  coupled  transport  phenomena  and 
their  repercussions  on  the  performance  and  lifetime  of  Li  ion  bat¬ 
teries. 
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Appendix  A.  List  of  symbols 


B  magnetic  induction 

D  diffusion  constant 

De>s  diffusion  constant  in  electrolyte  or  solid 

D  electric  displacement 

E  electric  field 

F  Faraday  constant 

H  magnetic  field 

M0,  M±  molar  mass  of  solvent,  Li  ion  and  counterions 
N±  flux  of  ions  or  counter  ions  (mol  m-2  s-1 ) 

Q.  renormalized  heat  flux  (see  Eq.  (17)) 

R  entropy  production  ( W  m-3 ) 

R  universal  gas  constant  (J  K-1  mol-1 ) 

T  temperature 

U0  open  circuit  potential 

c  Li  ion  concentration 

c±  concentration  of  Li  ions  and  counter  ions  (mol  m-3 ) 

ce,s  Li  ion  concentration  in  electrolyte  or  solid 

cp  specific  heat  per  unit  mass  at  constant  pressure 

e  energy  density  (J  m-3 ) 

ise  Butler  Volmer  flux 

z0  scale  of  Butler  Volmer  flux 

j  electrical  current 

k  Butler  Volmer  rate  constant 

1<b  Boltzmann  constant 

kT  Soret-Dufour  coefficient 

q  heat  flux  (Wm-2) 

s  entropy  density  (J  m-3 ) 

z±  number  of  elementary  charges  per  ion  or  counter  ion 

aa,c  apparent  anodic  and  cathodic  transfer  coefficients 
P  Seebeck  coefficient  (V  K-1 ) 
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k  ionic  conductivity 

X  heat  conductivity  (W  K-1  m-1 ) 

lx  collective  chemical  potential  of  Li  ions  (see  Eq.  (6) 

/x0,  ii±  chemical  potential  of  neutral,  positive  and  negative  phase 
/ it  Thompson  coefficient  (V  K-1 ) 

(p  electrochemical  potential  of  electrolyte 

p  mass  density  (kg  rrr3 ) 

pd  charge  density  (Cm-3 ) 

o  electronic  conductivity 

O  electrical  potential 

n  Peltier  coefficient  (V) 
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